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Abstract 

This paper presents a set of tools to compute topological information of sim- 
plicial complexes, tools that are applicable to extract topological informa- 
tion from digital pictures. A simplicial complex is encoded in a (non-unique) 
algebraic-topological format called AM-model. An AM-model for a given ob- 
ject K is determined by a concrete chain homotopy and it provides, in par- 
ticular, integer (co)homology generators of K and representative (co)cycles 
of these generators. An algorithm for computing an AM-model and the co- 
homological invariant HBl (derived from the rank of the cohomology ring) 
with integer coefficients for a finite simplicial complex in any dimension is 
designed here, extending the work done in jsl in which the ground ring was 
a field. A concept of generators which are "nicely" representative is also 
presented. Moreover, we extend the definition of AM-models to 3D binary 
digital images and we design algorithms to update the AM-model information 
after voxel set operations (union, intersection, difference and inverse). 
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1. Introduction 



The problem of adapting topology based metliods to discrete data is an 
active field of research. The set of algorithmic tools for computing topological 
properties in the setting of digital imagery is relatively small. Mainly, Betti 
numbers, Euler characteristic, skeletonization and local topological charac- 
terization. Recently, homology groups have also been computed in the digital 
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In order to enlarge this set, an algebraic-topological representation for 
a given geometric object is developed in ^ in terms of a chain homotopy 
equivalence (working with coefficients in a field). These models allows to 
compute the (co)homology groups, representative (co)cycles of the genera- 
tors, the cup product on cohomology as well as the cohomological number 
HBl. Our peculiar approach comes from the following main sources: (a) 
Eilenberg-Mac Lane work on homology of simplicial sets [sj; (b) Effective 



Homology Theory [23 



In this paper, we extend the results given in to the integer domain. 
Algebraic-topological representations of geometric objects are expressed here 
only in terms of particular chain homotopies (called AM- models). Compared 
to previous works, these models also encode torsion groups. Moreover, all the 
algorithms for computing integer homology based on the matrix reduction 
method to Smith normal form (for example [isl, [l|, 0, 19|) can be translated 



into our setting without extra computational cost. Finally, we successfully 
apply our computational algebraic-topological approach to 3D binary digital 
images. Moreover, AM-models for 3D binary digital images can be reused 
after voxel-set operations (union, intersection, difference and inverse). 

2. Algebraic-topological Models for Integer Homology Computa- 
tion 



In this section, we deal with the concept of AM-model (first given in 
and 1^) for the computation of all the integer homological information of a 



simplicial complex. In fact, we redefine this concept in simpler terms, what 
allows to store the same information in less space. 

We first introduce some basic algebraic-topological notions needed through- 



out the paper, which are extracted, mainly, from [18|. Let Z be the ground 
ring. Considering an ordering on a vertex set V, a simplex with q+1 affinely 
independent vertices Vq < ■ ■ ■ < Vq of V is the convex hull of these points, 
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denoted by cr'^ = {vq, . . . ,Vq). The dimension of a"^ is q. li i < q, an i-face of 
a'^ is an z-simplex whose vertices belong to the set {vq, . . . ,Vq}. A simplicial 
complex K is a collection of simplices such that every face of a simplex of 
K is also a simplex of K and the intersection of any two simplices of K is 
either a face of both or empty. The dimension of K is the one of the highest 
dimensional simplex in K. K^'^^ denotes the set of all the g-simplices of K. 
A chain complex C is a sequence 

... -^4^ ^ C,„i V ■■■ ^Co^O. 

of abelian groups Cq (called q-chain groups) and homomorphisms dq : Cq ^ 
Cg_i (called the differential of C in dimension g), indexed with the integers, 
such that dq-idq = for all q. In this paper, each group Cg is free of finite 
rank. A chain complex C can be encoded as a couple {C,d), where: (1) 
C = {Cq} and, for each q, Cq is a base of Cq] (2) d = {dq} and, for each q, dq 
is the differential of C in dimension q with respect to the bases Cq and Cq_i. 

Such an algebraic structure can be associated to a given simplicial com- 
plex K, in the following way: a q-chain is a finite sum ^ Xicrl where 
Aj e Z and af G K^'^\ The g-chains form the qth chain group of de- 
noted by Cq{K) (with g > 0). The boundary of a g-simplex = {vq, . . . ,Vq) 
is the (g — l)-chain dq{a'^) = Xli'=o(~-'-)'(^05 ■ ■ ■ ■ ■ ■ , Vq), where Vi means 
that Vi is omitted. By linearity, dq can be extended to g-chains. Then, the 
chain complex C{K) is the collection of chain groups Cq{K) connected by the 
boundary operators dq. 

Given a chain complex C = {C,d), a g-chain a'' G Cq is called a q- cycle 
if dq{a'') = 0. If a' = c/q+i(6«+^) for some G Cq+i then is called a 
q-boundary. Denote the groups of g-cycles and g-boundaries by Zq and Bq 
respectively. Define the qth homology group to be the quotient group Zq/Bq, 
denoted by l-lqiC). We say that two g-cycles a'^ and ¥ are homologous if 
there exists a (g + l)-chain c^"*"^ such that a'^ = ¥ + dq^i{c'^^^). For each g, 
the integer gth homology group l-iqiC) is a finitely generated abelian group 
isomorphic to Fq®Tg, where = Z © ■ ■ ■ © Z and = Z/tf © ■ ■ ■ © Z/f^^^ 
are the /ree subgroup and the torsion subgroup of 'Hq{C), respectively. The 
numbers t^, . . . , which satisfy that t\>2 and ^11^2! ' ' ' l^m, called the 
torsion coefficients of 'Hq{C). The rank of Fg, denoted by /3q, is called the gth 
Betti number of C. Intuitively, /3q is the number of components of connected 
pieces, (3i is the number of independent "holes" and /32 is the number of 
"cavities". For all g, there exists a finite number of elements of 'Hq(C) from 
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which we can deduce all 'Hq{C) elements. These elements are called homology 
generators of dimension q. We say that a'' is a representative q-cycle of a 
homology generator a'^ of dimension g if = + -Bg. We denote = [a'^]. 
Finally, the homology of a simplicial complex K is defined as the homology 
of C(i^). 

As far as the homology computation of a chain complex is concerned, 
the classical algorithm for computing homology with coefficients in Z is 



the integer reduction algorithm [18|]. Given a chain complex C = {C,d), 
this algorithm consists in reducing the matrix of the boundary operator in 
each dimension q, to its Smith Normal Form (SNF), relative to some bases 
{al, . . . o-mq} of {el~^, . . . , e^^^} of Cg_i such that for some tg, iq, Sq_i 

where 1 < Sq^i < mg„i and 1 < tg < ig < min(mg, 

• 9g_i(e^~^) = for 1 < ? < and (9g_i(e^~^) ^ for Sq_i < i < niq^i. 

• dg{a'^i) = for 1 < i < tg] 

• dg{a1) = Xyr\ where G Z and > 2 for < z < ig-, 

• and dq{al) = for £q < « < m^. 
In this case, 

= Z© ©Z and Tg_i = Z/A?^+i © ■ ■ ■ © Z/A^^ . 

Moreover, . . . , e^~_\} and {ejj!^;^, . . . , e^"^} are sets of representative 

cycles of the generators of and Tq_i, respectively. 

In [li, an algorithm improving the efficiency of this classical integer 



reduction algorithm is described. Their technique is mainly based on the 
results of 0] , in which a matrix reduction to integer SNF is determined in an 
efficient way. We can take advantage of these improvements in the algorithms 
described here without additional computational cost. 



A chain contraction 17| (/, 0) of a chain complex C = {C, d) to a chain 
complex C = {C',d') is a set of three homomorphisms / = {fg : Cg Cg}, 
g = {gg : Cg — t- Cg} and = {(pg : Cg Cg+i} such that for each q: 

• fq-idq = d'gfg and dggg = gg-id'g. 

• fggg is the identity map of Cg] 
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• (pq-idg + dg^i(f)g = 1(1 g " Qgfg, llidX Is, (f) IS &. chaiii homotopy of the 
identity map id = {idg} of C to gf. 

In this case, C has fewer or the same number of generators than C while C 



and C have isomorphic homology groups [l8|, p. 73]. 



A translation of the integer reduction algorithm in terms of chain con- 



tractions has been made in [8|, |10| • In those papers an AM-model for a given 
simplicial complex K was defined as a chain contraction of C{K) to a chain 
complex Ai such that the SNF, Ag, of the matrix of the differential of Ai 
in each dimension q satisfies that any non-null entry of Ag is greater than 1. 
Here we go further and define an AM-model for a given simplicial complex 
only in terms of the chain homotopy 0. We prove that it is possible to recover 
the other homomorphisms / and g and the chain complex A4. Moreover, our 
strategy outperforms the previous algorithms for computing integer homol- 
ogy in several points: (1) cohomological features can be computed; (2) we 
can control the topological changes after addition or deletion of simplices. 

Definition 2.1. An AM-model for a simplicial complex K is a couple {C, 0) 
such that: 

• C = {Cq} and, for each q, Cg is a base of Cq{K); 

• (f) = {(pq} and, for each q, (pg : Cq{K) — )■ Cg^i{K) is a homomorphism 
satisfying that (pq+i4>q = and = (pg; 

• the chain complex M. = {M,d) (such that Mg = Im iTg, Wq = idq — 
(pq-idg — dq+i(f)g and dg = dg\j^^^), satisfies, for each q, that any non- 
null entry of the SNF of dg is greater than 1. 

Theorem 2.1. Given an AM-model {C,(j)) for a simplicial complex K, we 
can define a chain contraction {f,g,(j)) of C{K) to a chain complex Ai such 
that any non-null entry of the SNF of the matrix of the differential of Ai, 
for each q, is greater than 1. In particular, if the homology of K is free then 
Ai is isomorphic to 'H{K). 

Proof. We only have to define / as {-Kg} where, for each g, tt^ = idg — 
(pq-idg — dq+i(f)q, aud g as the inclusion. Let us see that {f,g,(f)) is a chain 
contraction. It is clear that, for each q, gqfq = idg — (pq^idq — dg+i(f)g. On the 
other hand, let a"^ G Aiq. There exists b'^ G Cq{K) such that a'^ = i^q{b'^) = 

fe" - (pq-ldqib'^) - dq+^<Pq{¥). TheU, /,^7,(a'?) = 7r,(a'?) = ¥- (Pq^ldqib") - 
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dq+l(j)g(f)g_ldg{b'') + dg + l(l)gdg+l(f)g{bl) = fl^ . SeCOnCl, fg^lOg = g - 8 gCj) g^lO g = 

dgfg. Finally, if a'' G Mg then dg{a'^) E Mg. Therefore Qg-idg = dgQg. □ 
Given a simplicial complex i^', it is possible to define different AM-models 
for K since the chain homotopy (f) and the chain complex M. are not unique. 
In the following lemma, we show how to compute a new AM-model from a 
previous one. 

Lemma 2.1. Let [C, 0) be an AM-model for a simplicial complex K and let 
h'^eCgiK). 

• Ifhi e Mg then (t)g{h'>) = 0; and if dg{hi) = then 7rg{h'^) = /i"?. 

• Ifx'^e Cg{K) and h'^ e Mg such that dg{x'i) = and nq{x'^) = h'^ then 
[x'^] = [h'^] and we can define a new AM-model {C,(j)') as (f)'g{x'^) '■ = 
(f)g{h'^), := 4'q{x'^) and (f)' := (j) for the rest. In this case, vr^(a;^) = 
{idg - <f)'g^,dg - dg+i(j)'g){x'i) = and VT^ ( /l^ ) = x'^. 

Algorithm 2.2. Computing an AM-model for a Finite Simplicial Complex. 

Input: A simplicial complex K of dimension n. 
For q = 1 to q = n do 

reduce the matrix of dg to its SNF relative to some bases 
{a?,...,a^J of Cg{K) and {el', . . . , etJ^_J of such that 

for some where 1 < Sg_i < rrig^i): 

dg-i{el') =0, for 1 <i < Sg^i 
and dg-ile^') ^ 0, for Sg^i < i < rUg^i; 
for some tg,lg where ^ < tg < Ig < min(mq, : 
dqial) = elK for l<i<tg; 
dg{a1) = XlelK XjeZ, \1>2 for tg<i< ig-, 
and dg{al) = for ig < i < nig . 
Define Cg^^ := {el', . . . ,6^1^} , := {a?, . . . , a^J , 
4>g-i{e1 ') := a] for l<i<tg, 
0g_i(e^"^) := for tg<i< Sg^i 
and 0g(a^) := for 1 <i <tg. 
Output: The couple (C, 0) . 

Notice that, for q = k, 4>k{o,'^) is defined for 1 < z < tk- After, (j)k is again 
defined when q = k -\- 1, but not for the previous elements , since ^ Im 
dk+i, for 1 < z < tk. 
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Theorem 2.2. Let K be a finite simplicial complex of dimension n. The 
output of Algorithm \2.^ {C,(j)), defines an AM-model for K. Moreover, 
integer homology generators and representative cycles of these generators can 
he directly obtained from Ai. 

Proof. First of all, let us see that for each g, if x'' G Cg, then 7Tq{x'^) = 
or Hg{x'^) = x'^, where vr^ = idg — (pg-idg — dg+icpg. Since the matrix of the 
differential in dimension q with respect to the base Cg coincides with its SNF, 
then for each G C^: 

• If dg{x'^) = y^~^ for some y'^^^ G Cg_i, then (j)g^i{y'^~^) = x'' and 

= 0. Therefore T^g{x^) = 0. 

• If dg{x'') = Ay«-\ for some y^"^ G Cg^i and A G Z, A > 2, then 
(f)q-i{y'^~^) = and 0q(x^) = 0. Therefore, Hq^x'^) = x^. 

• If dg{x'^) = and there exist z"'^^ G Cg+i such that dg+i{z'^~^^) = x'^ 
then = and iTgix") = 0. 

• If dg{x'^) = and there exist z^"*"^ G C^+i and A G Z, A > 2, such that 
dg+i{z'^+^) = Xx", then (f)g{x'^) = and iTgix") = x". 

• If dg{x'^) = and there is neither A G Z, A 7^ 0, nor z'^'^^ G C^+i, such 
that = Ax", then (pg^x'^) = and 7ig{x'^) = x'^. 

Therefore, a base of M. in each dimension q is the set Mg = {x" : x'' G 
and vrg(x") = x"}. Now, for each q, let {x^, . . . ,x^_^} be the elements of Mg 
and . . . , y^i^} the elements of Mg_i. For some 1 < Sg^i < nig^i, 

dg^i{y1~'^) = for 1 < z < Sg_i and dg^i{y^'~^) ^ for Sg^i < i < mg_i. For 
some ig where I < ig < min(mq, Sg_i), 

• dgix'^i) = Xlyr\ where A^ G Z and A^ > 2 for 1 < z < ig-, 

• dg{a1) = for £q < i < mg. 
In this case, 

= Z© ©Z and T^^i = Z/\\ © ■ ■ ■ © Z/A^^ . 

Moreover, {y'^'^^, . . . , and {i/i""*^, . . . , yl~^} are sets of representative 

cycles of the generators of Fg_i and Tg_i, respectively. □ 
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Figure 1: The Klein bottle and a triangulation of it. 



Example 2.3. Consider the simplicial complex K in Figure 1 whose un- 
derlying space is the Klein bottle \la . p. 283]. Applying Algorithm \2.2\ we 
obtain that the vertex (a) is an element of Cq. The rest of the elements of 
Co are the boundaries of the 1-simplices marked in blue in Figure 1. These 
1-simplices are also elements of Ci. Denote by x one of these 1-simplices. 
The cycles ai := {a,d) + {d,e) — (a, e) and a2 := {a,b) + (6, c) — (a, c) are 
also elements of Ci . The rest of the elements of Ci are the boundaries of all 
the 2-simplices except for {f,g,h). These 2-simplices belong to C2. Denote 
by y one of these 2-simplices. The 2-chain consisting in the sum of all the 
triangles in K, 7 := -{a,b,f) - {b,c,f) + {a,c,g) - {a,e,g) + {e,g,i) - 
(e, d, i) + (c, d, i) — {a, c, d) + {b, c, i) + (a, b, h) — (a, e, h) — (e, d, f) + (a, d, f) + 
(c, /, g) — (/, g, h) — {h, g, i) — {b, h, i) , is also an element of C2. The image of 
d, and n on C and the SNF of the matrix of the differential in dimension 
1 and 2 are given below: 



c 


(a) 


di{x) 


X 


ai 


a2 


d2{y) 


y 


7 


d 








di{x) 











d2{y) 


2ai 








X 











y 








TT 


(a) 










a2 








7 





X 


d2{y) 




^2 


di{x) 


1 











(a) 















d2 


y 


7 


d2{y) 


1 





ai 





2 


02 









{01,02} and M2 = {7}, d\Mo{{a)) = 0, 
,(7) = 2ai. We obtain that 1-Lq{K) ~ 



Therefore, Mq = {(a)}, Mi 

= 0, 9|a4i(«2) = and d\M2{l) = 2ai 
Z, T-Li{K) ~ Z © Z/2 and representative cycles of the homology generators 
are (a) for 'Hq{K), a2 for the free part of'Hi{K) and ai for the torsion part. 
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3. Extracting Integer Cohomology Information from AM-models 

In this section, we extend tlie work done in jof (witli coefficients in a 
field) for computing cofiomology (tfie dual notion of liomology) over the 
integer domain. Working with coefficients in a field, homology groups are 
free and isomorphic to cohomology groups. Nevertheless, integer homology 
and cohomology can have a torsion part and they are not isomorphic, in 
general. The cohomology groups have an additional multiplicative structure, 
the cup product, from which we can derive finer invariants than homology. 

Let C = {C,d) be a chain complex. The cochain complex C* = {C*,6) 
in each dimension q is the group of q-cochains with coefficients in Z, = 
{c : — )■ Z such that c is a homomorphism}. If Cq = {oi, . . . , amq} is a 
base of Cq then a base of is = {al, . . . , a^q}' "where a* : — j- Z is 
given by ^^(aj) = 1 and a*{aj) = for 1 < z,j < rriq and j ^ i. For all 
q, the differential dq^i on C^+i induces the codifferential 5'^ : ^ C^^^ via 
5'^{c) = cdq+i, so that S'^ increases the dimension by one. Define Z*^ to be the 
kernel of S'^ and B'^^^ to be its image. These groups are called the group of 
q-cocycles and q + 1-coboundaries, respectively. Define the gth cohomology 
group, 'H'^{C) = Z'^/B'^ for g > 0. For each q, the integer gth cohomology 
group WiC) is a finitely generated abelian group isomorphic to F'' (B T'^, 
where F'^ and T'^ are the free subgroup and the torsion subgroup of 'H'^{C), 
respectively. The rank of F'^ coincides with the rank of Fq for each q. 

Theorem 3.1. Let {C,(f)) be the output of Algorithm \2.^ for a given finite 
simplicial complex K. The integer cohomology of K and representative cocy- 
cles of integer cohomology generators can be directly obtained from M. = im 
IT, where -k = id — (pd — dcj). 

Proof. A base of M. is the set M = {x : x & C and ii{x) = x}. Now, for each 
q, let {xi, . . . , Xmq} be the elements of Mq and {yi, . . . , ymq-A the elements 
of Mg_i. For some 1 < Sg-i < mq_i, dq-i{yi) = for 1 < < Sq-i and 
dq-iiiji) 7^ for Sq-i < i < niq^i. For some ig where I < £q < min(mq, Sq-i), 
dq{xi) = Xiyi, where Aj G Z and Aj > 2 for 1 < i < iq] and dq{xi) = for 
iq < i < rriq. In this case, 

F^-i = Z© ©Z and = Z/Ai © ■ ■ ■ © Z/\ . 

Moreover, • • • iy*s^_^} and {x\, . . . ,x}^} are sets of representative co- 

cycles of the generators of F'^~^ and T'^, respectively. □ 
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Example 3.1. Consider the AM-model {C,(f)) obtained in Example \2.3\ for 
the simplicial complex K whose underlying space is the Klein bottle. Starting 
from the chain complex M. whose base is {(a), ai, 02, 7} and differential d\j^, 
we construct in a straightforward way the cochain complex Ai* whose base 
is {(a)*, aj, ttj) 7*} o-nd codifferential 5 given by: 6o{{a)*) = {a)*d\Mi ~ ^' 
5,(al) = ald\^^ = 0, 5,{a*,) = ald^^ = 2^ , Hi*) = 1*9^, = 0- 

Therefore we obtain that n"{K) ~ Z, V}{K) ~ Z and n^{K) ~ Z/2; 
and the representative cocycles are: {a)* in dimension 0, al in dimension 1 
and 7* in dimension 2. 

The cochain complex C*{K) is a ring with the cup product C^{K) x 
C^iK) -> CP+'^{K) given by: 

This product induces an operation T-L^i^K) x T-U^i^K) — )■ 'W^''{K), via 
[c] [d] = [c ^ c'], that is bilinear, associative, commutative up to a sign, 
independent of the ordering of the vertices of K and homotopy-type invariant 



18|, p. 289]. 

Working with coefficients in Z/2, a new cohomology invariant called HBl 
is obtained in [9|. The idea is to put into a matrix form the multiplication 
table of the cup product of cohomology generators of dimension 1. The 
following algorithm compute HBl working with integer coefficients. Assume 
that an AM-model for /C, (C, 0), is computed using Algorithm 12. 2[ Then, 
for each g, Mq = {x : x E Cq and vrg(x) = x} is a base of A^g = im vr^, 
where tt, = idq - (j)qdq - dq(j)q. Let = {x* : x e Mg} where x* : Mq Z 
is such that x*{x) = 1 and x*{z) = for 2; G Mg and z ^ x. Suppose that 
there are r, elements in and Sj cocycles in M*, i = 1, 2. Then the following 
algorithm computes the cohomological invariant HBl, working with integer 
coefficients, in C(rJs^r2S2). 

Algorithm 3.2. Algorithm for computing the cohomological invariant HBl. 

Input: An AM-model (C, 0) for a simplicial complex K computed 
using Algorithm 12 . 21 . 

Let M* = {x* : x e C and 7r(x) = x} . 

Let {al, . . . ,a*} and {ll, ■ ■ ■ ,1^} ^® ^^^^ °^ 1 2-cocycles 
in M* , respectively. 
For 1 = 1 to p do 
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Figure 2: The simplicial complex T, representative cocyclcs of the generators of 1H}{T) 
and the muhiphcation table of the cup product. 



For j = i to p do 

For A; = 1 to m do 

:= (a* -a*)(7fc). 
HBl := the rank of the 2D matrix of integers {b[[ij)^k)) ■ 
Output: The integer HBl. 



The implementation of the algorithm described above has already been 
made 
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We have tested it on several 3D objects. We give here an example 
of the computation of the cohomology, representative cocycles of cohomology 
generators and the invariant HBl. 



Example 3.3. Consider the simplicial complex T whose underlying space 
is showed in Figure 2 (on the left). It consists in 11847 simplices. The 
running time for computing an AM-model for T and the homology ofT using 
a Pentium 4, 3.2 GHz, 1Gb RAM is 2 seconds. We obtain that (Sq = 1, f3i = 4 
and = 3. The running time for computing the cup product is 1.5 seconds. 
In Figure 2 (on the center), the 1 and 2-simplices, on which the representative 
cocycles are non-null, are drawn. The table on the right of Figure 2 shows 
the results of the cup product of any two cohomology generators of dimension 
1. Finally, HBl = 2. 



4. AM-models for 3D Digital Images 



An important issue in Digital Volume Processing is to design efficient al- 
gorithms for analysis and processing in grids such as the face-centered cubic 
(fee) and the body-centered cubic (bcc) grids since it is very easy to 
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obtain data structures for them. The bcc and the fee grid are the general- 
izations to 3D of the two-dimensional hexagonal grid. In the bcc grid, the 
voxels consist of truncated octahedra, and in the fee grid, the voxels consist 
of rhombic dodecahedra. They are better approximations of Euclidean balls 
than the cube. 

A 3D digital image X is encoded as a tuple (V, /, b, w), where V is the set 
of grid points in a 3D grid, / is the set of black points and b (resp. w) deter- 
mines the neighborhood relations between black points (resp. white points) 
in the grid. A bcc grid is equivalent to a grid V in which the grid points are 
those (xi,X2,X3) G such that xi = X2 = X3 (mod 2) (see [l5|). The only 
Voronoi adjacency relation on V is the 14-adjacency. Using this adjacency, it 
is straightforward to associate to a digital image X = (V, /, 14, 14) a unique 
simplicial complex K{I) (up to isomorphism) with the same topological in- 
formation as /. It is called the simplicial representation of I. The i-simplices 
of K{I) {i G {0,1, 2, 3}) are constituted by the different sets of i mutually 14- 
neighbor points in /. Since an isomorphism of digital images is equivalent to 
a simplicial homeomorphism of the corresponding simplicial representations, 
we define the (co)homology of / as the (co)homology of K{I). Moreover, we 
define an AM-model for X as an AM-model for its simplicial representation 
K{I), {Ci,(j)i), where Ci is a base for C{K{I)). Therefore, the simplicial 
complexes considered in this section are embedded in R^, then their homol- 
ogy groups vanish for dimensions greater than 3 and they are torsion-free for 
dimensions 0, 1 and 2 (see 0, ch.lO]). Moreover, the value of all the possible 
non-null entries of the SNF of the matrix of the differential of C{K) (where 
K is a simplicial complex embedded in R^) in each dimension must be 1; 
and an AM-model (C, (p) for K satisfies that the chain complex = im 
71 = id — (f)d — d(f) is isomorphic to the homology of K. 



All the algorithms explained in this paper have been implemented [11 
using as a grid the set of points with integer coordinates in the Euclidean 
3-space and the 14-adjacency by which the neighbors of a grid point 
(black or white) with integer coordinates (xi,X2,X3) are: (xi ± 1,X2,X3), 

(Xi, X2 ± 1, X3), (Xi, X2, X3±l), (Xi + 1, X2 - 1, X3), (Xi - 1, X2 + 1, X3), (Xi + 
1,X2,X3 - 1), (Xi - 1,X2,X3 + 1), (Xi,X2 + 1,X3 - 1), (Xi,X2 - 1, X3 + 1), 

(xi -|- 1, X2 -l- 1, X3 — 1), (xi — 1, X2 — 1, X3 + 1). The digital spaces {Z^, 14, 14) 
and (V, 14, 14) are isomorphic: a grid point (xi, X2, X3) of (Z^, 14, 14) can be 
associated to the point (xi + X2 + 2x3, — xi + X2, — xi — X2) of (V, 14, 14). 

In order to compute an AM-model for a digital image X = {Z^, 1, 14, 14), 
we take advantage of the particular structure of K{I) in the way that we 
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Figure 3: The 14-neighbors of a point p (on the left) and p G V (on the right). 



consider as the input of Algorithm 12.21 the following special initial base: 

• In dimension 0, it consists in the set of the vertices of K{I), except for 
the vertices {{xi,X2, x-s)) such that X3 is odd, which are replaced by the 
boundary of the 1-chains a = ((xi, X2, X3 — 1), (xi, X2, X3)), if a G K{I). 

• In dimension 1, it consists in the set of all the edges in K{I), except for 
the edges of the form ((xi, X2, X3), (xi, X2 + I, a^s)), which are replaced by 
the boundary of the 2-chains b = ((xi, X2, X3), (xi, X2, X3 + 1), (xi, X2 + 
1,X3)) if 6e J^/). 

• In dimension 2, it consists in the set of all the triangles in K{I), except 
for the triangles of the form ((xi — 1,X2 + 1,2:3), (xi,X2,X3), (xi,X2 + 
1,X3)), which are replaced by the boundary of the 3-chains c = ((xi — 

1,X2,X3 + 1), (xi-l,X2 + l,X3), (xi,X2,X3), (xi, X2 + 1, X3)), if C G K{I). 

• In dimension 3, it consists in the set of all the tetrahedra in K{I). 

We then reduce the matrix of dq relative to this base to its SNF and it holds 
that we do not have to modify the rows and columns corresponding to the 
chains a, 6, c, (9(a), d{h) and d{c). 

In the following table we present the running time for computing AM- 
models for the 3D digital images showed in Figure 4 using the program 



developed in [H 



Image X 


Number of points of / 


Time for computing 


/3o 


/3i 


/32 


h 


26308 


50 seconds 


2 


9 


3 


h 


31012 


38 seconds 


138 


419 


13 


h 


18842 


27 seconds 


1 


277 


5 
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Figure 4: The 3D digital images /i, I2 and I3. 

4-1 ■ Computing "Good" Representative Cycles of Homology Generators 

In [gI , algorithms for obtaining "optimal" generators of the first homology 
group are developed using Dijkstra's short est p ath algorithm for any oriented 



2-manifolds. The algorithms presented in |20[ produce "nice" representative 
cycles of homology generators that always fit on the boundaries of the image. 

We sketch here some techniques for drawing "good" representative cycles 
of homology generators in the context of digital volumes. Let (C^, (pi) be an 
AM- model for a 3D digital image X = (V,/, 14, 14). We denote by dX = 
(V, dl, 14, 14) the digital image such that dl is the set of points of / with at 
least one 14-neighbor white point of X. We say that M is a set of "good" 
representative cycles of homology generators of I if each chain of M satisfies 
that: 



• it is a cycle; 

• it belongs to C{K{dI))- 

• in dimension 0, it is a vertex; 

• in dimension 1, it is an elementary cycle (that is, it is connected, each 
vertex is shared by exactly two edges and two consecutive edges can 
not belong to the same triangle in K{I)); 

• in dimension 2, it is an elementary cavity (that is, it is a connected 
2-cycle with exactly one white connected component inside and three 
triangles can not belong to the same tetrahedra in K{I)). 

First of all, suppose we have an AM-model for dX, {Cei, (j^ai), and X, (Cj, 0j). 
Let iTgj = id — d(f)gi — (pgjd and ttj = id — d(f)j — (pid. Let us denote by 
{«!, . . . , an} the elements of Mg^ which is a base of JUgi = im Hgi. Suppose 
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Figure 5: Representative cycles of the homology generators of the images /i, I2 and 1^. 

that if e Mqj then d{h) = (we can obtain this if we compute the AM- 
model for dl using Algorithm 12. 2p . The cycles of Mgi are representative 
cycles of the homology generators of dl. Decompose and replace each 0- 
cycle in Mgj by its constitutive vertices, each 1-cycle in Mgj by its elementary 
cycles and each 2-cycle in Mgj by its elementary cavities. Note that all the 
homology generators of / are homology generators of dl, therefore a set of 
representative cycles of the homology generators of / is a subset of Mgj. Let 
M'j be the set {a : a G Mgj and ni{a) = a}. Obtain the new AM- model 
(Cj, applying Lemma [211 to the elements of the set M^. Then, (Cj, 0^) is 
a new AM-model for X; and Mj is a base of A^j = im tTj and a set of good 
representative cycles of homology generators of /. 

4-2. AM-models after Adding or Deleting a Voxel 

Now, we study the problem of topologically controlling a digital image 
using AM-models when it suffers local changes. More concretely, we show 
how to compute an AM-model for a digital image when a voxel is added or 
deleted using the AM-model computed before. Adding or deleting a voxel v 
of I means to change the color of a point f in X and it consists of adding or 
deleting a set of simplices having f as a vertex. Since we work with simplicial 
complexes representing 3D digital images considering the 14-adjacency, the 
maximum number of simplices having as a vertex is 74. The key idea of 
both algorithms is that when a g-simplex a is added to or deleted from an 
AM-model for a simplicial complex K, we only have to put into the SNF the 
column of the matrix of dq relative to a to obtain the new AM-model. 
AM-models after Adding a Voxel. Given a digital image X = (V, /, 14, 14), 
we add to / a point v ^ I to obtain a new digital image denoted by 
jy^' = (V,/ U {f},14, 14). Therefore, the addition of a point f to J con- 
sists in the addition to K{I) of all the simplices of K{I U {v}) having v as 
a vertex. In each step of the process, one simplex is added. Given an AM- 



IS 



model for a digital image X = (V, /, 14, 14) such that / has m points, the 
following algorithm computes an AM- model for X^*', with integer coefficients, 
in 0{rn?) 

Algorithm 4.1. An Incremental Algorithm for Computing an AM-model 
for a 3D Digital Image X. 

Input: A digital image (V, /, 14, 14) , a point v ^ I and an AM-model 
{Ci,(j),) for X such that, in each dimension q, the 
matrix of dq with respect to C, coincides with its SNF. 

Let {o"i, . . . , (T„} (n < 74) be the ordered-by-increasing-dimension 

set of all the simplices of K{I {v}) having f as a vertex. 
:= and := 0, . 

For i = 1 to i = n do: 

Let q be the dimension of at, let C^^ = {ai, . . . , a^} and 

Cg!^i = {ci, . . . , e^} such that dg{aj) = Cj for 1 < j < t , 

dg{aj) = for t < j <r and dg{ai) = J2i=iP(^^£ some G Z . 

a := (Jj - XlLi Pe^f^ «r-, a} ■ 

If = for £ > t then (f)^'"{a) := 0. 

Else obtain the SNF of the matrix of dq 

relative to some base C^l^ := {ei, . . . , ef, e^.,.]^, . . . , e'^} then 
0Ut.(a) .- 0^ 0U''(e;+i) := a and 0^''(e;.) := for the rest. 
Output: An AM-model (C^'',0^'') for . 

AM-models after Deleting a Voxel. Given a digital image X = (V, /, 14, 14), 
we delete from / a point f G / to obtain a new digital image denoted by 
X\*' = (V, /\ {t"}, 14, 14). Therefore, the deletion of a point v from I consists 
in the deletion from K{I) of all the simplices of K{I) having as a vertex. 
In each step of the process, one simplex is deleted. Given an AM-model 
for a digital image X = (V, /, 14, 14) such that / has m points, the following 
algorithm computes an AM-model for X^", with integer coefficients, in 0{m?) 

Algorithm 4.2. A Decremental Algorithm for Computing an AM-model for 
a 3D Digital Image X. 

Input: A digital image (V, /, 14, 14) , a point vel and an AM-model 
{Ci,(f)j) for X such that in each dimension q, 
the matrix of dq with respect to Cj coincides with its SNF. 



16 



Let {o"i, . . . , o"„} (n < 74) be the ordered-by-decreasing-dimension 

set of all the simplices of K{I) having f as a vertex. 

K := K{I) , C^" := C, and := 0, . 

For 2 = 1 to i = n do 

Let q be the dimension of CTj, let Cg" = {«!,..., a^} and 

C^_^ = {ei, . . . , e^} such that dg{aj) = ej for 1 < j < t, dq{aj) = 

for t < j <r and 9g(cr,j) = X]|=i P^^^ some G Z. 

Let k be the smallest index such that {oi, . . . , a^, . . . , a^.} is 

a base of Cq{K \ {cXi}) then 

C}'" := {ai, . . . , 4, . . . , a J and 0(efc) := . 
K:=K\{a,}. 

Output: An AM-model for X^^ 

Observe that if an AM-model (C, (p) for X is computed using Algorithm 
\2.2\ l4.1l or I4l2] then in each dimension q, it satisfies that the matrix of dq with 
respect to C coincides with its SNF. 

4-3. AM-models for 3D Digital Images under Voxel-Set Operations 

In this subsection, we reuse the AM-model information for digital images 
under voxel-set operations (union, intersection, difference and inverse). Let 
X = (V, /, 14, 14) and J = (V, J, 14, 14) be two digital images, then XU J = 
(V, / U J, 14, 14), X n = (V, / n J, 14, 14), X \ J = (V, / \ J, 14, 14). In 
order to define the inverse of X, for each p = {xi,X2,Xs) G V, let Xp = 
max{|a;j| : i = 1,2,3}. Let Xj = ma.x{Xp : p G /}. Let Qj be the digital 
image (V, Gi, 14, 14) where Gi = {p : p eV and Xp < Xj + 1}. Then, the 
inverse of X, X, is Qi \X. We will not consider any of the trivial cases / = 0, 
J = 0, / nJ = 0, JCJ, or/ = {j9: pGV and Xp < r} for some r G Z. 

Let C = (V, L, 14, 14) be a digital image and F = {vi, . . . , Vm} C V such 
that F C L or FnL = 0. If F C L, denote by the image (V, L\F, 14, 14). 
On the other hand, if F fl L = 0, denote by the image (V, L U F, 14, 14). 
Let [Cl, (Pl) be an AM-model for C such that in each dimension q, the matrix 
of dg with respect to coincides with its SNF. Algorithm 14.31 is a common 
processing to the four voxel-set operations treated here. 

Algorithm 4.3. Common Processing. 

Input: The AM-model {Cl,4>l) for £ = (V, L, 14, 14) and the set of 
points F = {vi, . . . , Vm} such that F G L or FnL = 0. 
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U F C L then 

For i = 1 to i = m do 

apply Algorithm 14.21 to Vi and the AM-model {Cl,(()l). 
C,:=CY' and 0, := 
Else for i = 1 to i = m do 

apply Algorithm 14. II to Vi and the AM-model (C^, 0^) . 
:= C^^'' and 0^ := 0^^', . 
Cf := and 0f := 0^ 
Output: An AM-model (Cf,0f) for 

Let (Cj, 0j) and (Cj, 0j) be an AM-model for X and J', respectively, such 
that in each dimension q, the matrix of dq with respect to C^, for L = I, J, 
coincides with its SNF. The following algorithm computes an AM-model for 
lUj. 

Algorithm 4.4. Computing an AM-model for X U J' . 

Input: The AM-models (Cj,0j) for I and (Cj,0j) for J. 
Apply Algorithm HH] to (Ci,0£) and F := I n J for L = I,J. 
For each a G C := C7f U 

0(a):=0f(a) if aGCf; and 0(a) := 0^(a) if a G CJ . 
Apply Algorithm 1473] to F and the AM-model (C, 0) . 
C/uj := and 0^uj := 0"" ■ 

Output: An AM-model (Cjuj,0/uj) for I Li J. 

Algorithm 14. 31 is the essential step for computing an AM-model for XfljT, 
X \ and X. 

Algorithm 4.5. Computing an AM-model for X (1 J'. 

Input: The AM-model (Cj,0^) for X and the set F = I\J. 
Apply Algorithm 1473] to (C^,0j) and F. 
Output: An AM-model (Cf,0f) for I D J . 

An algorithm for computing an AM-model for X \ (resp. for X) is 
similar to the one above. The only difference is that the input is an AM- 
model (Cj, 0j) for X and the set F = / fl J (resp. an AM-model {Cqj, 4>Gi) 
for Qi and the set /). 
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5. Comments and Future Work 



The (non-unique) algebraic-topological representation of a given sim- 
phcial complex of any dimension showed here, allows to compute integer 
(co)homology, representative cycles of integer (co)homology generators, the 
cup product on cohomology with integer coefficients and a topological invari- 
ant derived from the integer cohomology ring. Moreover, we give a positive 
answer to the problem of reusing AM-models for determining homological 
information of new 3D binary digital images constructed from the previous 
ones using voxel-set operations. 

There is considerable scope for further research: 

• To extend our method to nD binary digital images in any grid using 



simplicial analogous techniques [IJ, ll3|, y 



To compute topological invariants from primary cohomology operations 
in the discrete setting of digital images. The work done in jsj seems to 
be a compulsory reading for advancing in this issue. 

To compute homotopy groups of nD binary digital images or n-G-maps 
191 using AM-models. The works [l6|, Q, El, H Q could help us m 



this task. 



To develop a discrete Morse theory 0] for digital images well-adapted 
to our method. The paper 22|] would be a good starting point. 



Potential applications of our method in Computer Vision and Digital 
Image Processing involving not only 3D object but also higher dimensional 
structures can be encountered in Medical Imaging and Object Modelling. 
Moreover, our method seems to be especially well-adapted to segmentation 
under topological constraints and elimination of small topological noise. 



References 

[1] Agoston M.K.: Algebraic Topology, a First Course. Marcel Dekker Ed., 
1976 

[2] Alexandroff P., Hopf H.: Topologie I. Springer, Berlin 1935 

[3] Ayala R., Dommguez E., Frances A.R., Quintero A.: Homotopy in Dig- 
ital Spaces. DGCIOO, LNCS 1953 (2000) 3-14 



19 



[4] Dumas J.G., Saunders B.D., Villad G.: On Efficient Sparse Integer Ma- 
trix SNF Computations. J. of Symbolic Computation 32 (2001) 71-99 

[5] Eilenberg S., MacLane S.: On the Groups H[7r,n), I and II. Annals of 
Math., 58 (1953), 55-106 and 60 (1954) 49-139 

[6] Erickson J., Whittlesey K.: Greedy Optimal Homotopy and Homology 
Generators. Proc. of the 16th Symp. on Discrete Alg. (2005) 1038-1046 

[7] Forman R.: A Discrete Morse Theory for Cell Complexes. Geometry, 
Topology and Physics for Raoul Bott, S.T. Yau, ed.. International Press, 
1995 

[8] Gonzalez-Di'az R., Real P.: Computation of Cohomology Operations 
on Finite Simplicial Complexes. Homology, Homotopy and App. 5 (2) 
(2003) 83-93 

[9] Gonzalez-Di'az R., Real P.: On the Cohomology of 3D Digital Images. 
Discrete Applied Math. 147 (2005) 245-263 

[10] Gonzalez-Diaz R., Medrano B., Real P., Sanchez-Pelaez J.: Reusing 
Integer Homology Information of Digital Images. DGCI06, LNCS 4245 
(2006) 199-210 

[11] Gonzalez-Di'az R., Sanchez-Pelaez J., Real P.: VOXELO Vi- 
sualizing Topology Computations on Binary Voxel Volumes, 
http: / /imagen-a.us.es 

[12] Herman, G.T.: Geometry of Digital Spaces. Birkhauser, Boston, 1998 

[13] Khalimsky E., Kopperman R., Meyer P.: Computer Graphics and Con- 
nected Topologies on Finite Ordered Sets. Top. and Appl. 36 (1990) 
1-17 

[14] Kong T.Y., Roscoe A.W., Rosenfeld A.: Concepts of Digital Topology. 
Topology and its App. 46 (1992) 219-262 

[15] Kovalevsky V.A.: Discrete Topology and Contour Definition. Pattern 
Recognition Letter 2 (1984) 281-288 

[16] Malgouyres R.: Computing the Fundamental Group in Digital Spaces. 
IJPRAI 15 (7) (2001) 1075-1088 



20 



MacLane S.: Homology. Classic in Math., Springer- Verlag, 1995 

Munkres J.R.: Elements of Algebraic Topology. Addison- Wesley Co. 
1984 

Peltier S., Alayrangues S., Fuchs L., Lachaud J.: Computation of Ho- 
mology Groups and Generators. Computers and Graphics 30 (2006) 
62-69 

Peltier S., Ion A., Haxhimusa Y., Kropatsch W.G., Damiand G.: Com- 
putingf Homology Group Generators of Images Using Irregular Graph 
Pyramids. GbRPR 2007, LNCS 4538 (2007) 283-294 

Real P.: Homological Perturbation Theory and Associativity. Homology, 
Homotopy and Applications 2 (2000) 51-88 

Skoldberg E.: Morse Theory from an Algebraic Viewpoint. Trans. Amer. 
Math. Soc. 358 (2006) 115-129. 

Sergeraert F.: Homologie Effective, I and II. C. R. Acad. Sci. Paris Ser. 
I Math. 304 (11) (1987) 279-282 and 304 (12) (1987) 319-321 



21 



